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We have constrained the charge-mass (e - m) phase space of millicharged particles through the simulation of 
the rotational evolution of neutron stars, where an extra slow-down effect due to the accretions of millicharged 
dark matter particles is considered. For a canonical neutron star of M = 1.4 Mq and = 10 km with typical 
magnetic field strength So = lO'^ G, we have shown an upper limit of millicharged particles, which is compatible 
with recently experimental and observational bounds. Meanwhile, we have also explored the influences on the 
s - m phase space of millicharged particles for different magnetic fields So and dark matter density in the 
vicinity of the neutron star. 
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I. INTRODUCTION 

Since Zwicky (1933) proposed the problem of the “missing 
mass”, the theoretical and experimental studies on dark mat¬ 
ter (DM) have attracted more and more attention. It is well- 
known that the total energy of the Universe contains about 
27% DM and about 5% baryon matter as well as about 68% 
dark energy according to the most recent cosmological re¬ 
sults based on Planck measurements of the cosmic microwave 
background (CMB) temperature and lensing-potential power 
spectra . In recent years, experimental physicists have con¬ 
centrated their attention on the direct and indirect detections 
of DM particles. Direct search experiments aim to detect indi¬ 
vidual interactions, i.e., DM particles scatter off target nuclei 
of detectors, such as DAMA S, CDMS II Ot], CoGeNT |01, 
CRESST-II lH, XENONIOO ll, LUX H, PandaX H and so 
on. Indirect detection experiments mainly probe annihilation 
or decay products (gamma-rays, neutrinos and charged cos¬ 
mic rays) of DM, such as AMS-02, PAMELA, Eermi-LAT, 
lACTs, IceCube, ANTARES, Super-K IKH. Erom a theo¬ 
retical point of view, several candidates of DM particles are 
supposed, for example, neutralinos in El, Majorana neutri¬ 
nos and technibaryons IH El . In this paper, we are inter¬ 
ested in another candidate of DM in the form of millicharged 
(MC) particles iEsUH. 

The electric charge of all the particles in the Standard 
Model appears to be an integer multiple of e/3, where e is 
the charge of the electron. However, MC particles have elec¬ 
tric charge e' - ee, where e is any real number and e < 1. 
They can be either bosons or fermions. MC particles were 
first proposed in order to solve the DM puzzle a long time 
ago lEx EM El . Hereafter millicharged dark matter particles 
(MCDM) have been studied extensively (see Ref. | E^ and 
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references therein). It is worth noticing that Huh et al. 
presented a possible explanation of the 511 KeV galactic y- 
ray due to positron productions from MCDM. MC particles 
have an obvious impact on the standard picture of the Uni¬ 
verse in many ways. Eirstly MC particles can significantly 
influence the expansion rate of the Universe and the baryon- 
to-photon ratio during the epoch of big bang nucleosynthesis 
(BBN) in the early Universe lEhll. Secondly, MC particles 
also may explain the creation of galactic magnetic fields at 
the cosmological epoch of the galaxy formation lEdll . Thirdly, 
the anisotropy power spectrum of the CMB is affected by MC 
particles in several respects [EHE^- 

At present MC particles are already constrained by experi¬ 
mental and observational data iIEHE^IztU^ . including the 
e-m (m is the mass of MC particle) phase space and the frac¬ 
tion of MC particles. These existing constraints on MC par¬ 
ticles mainly contain laboratory bounds (inclu ding the limits 
from experiments searching for MC particles cos¬ 

mological bounds (BBN, CMB anisotropy) and constraints 
from stellar evolution (globular clusters, white dwarfs, red gi¬ 
ants, and supernova 1987A). Eor very light MC particles, m < 
1 keV, the reactor experiments show a strict limit e < 10*^ 
iH. Eor MC particles that are lighter than electrons, m < m^, 
the best particle physics bound e < 3.4x 10“^ is given from the 
data that ortho-positronium decay to dark matter pairs j^ . 
Eor heavier MC particles, m > lUg, the bounds become weaker 
Eor m — 1 MeV, the bound is e < 4.1 x 10“^, and for 
m - 100 MeV, the bound is up to e < 5.8 x 10 If the MC 
particles become heavier than 100 MeV, the bound e ~ 10“^ 
is allowed, while for m> \ GeV it can be as large as e = 0.1. 

With the discovery of neutron star (NS) by Bell and Hewish 
IHt], the properties of NS have been studied by many theoret¬ 
ical and experimental physicists ll33 - [j^ . NS is a compact 
star, which is characterized by strong gravitational field, elec¬ 
tromagnetic field, extreme strong and weak interaction. Thus 
it is usually regarded as a “natural laboratory” with extreme 
physical condition. At present it seems that all currently ob¬ 
served pulsar periods mainly lie between milliseconds and a 
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few seconds. As it can be seen from the P-P (P and P are the 
spin period and its time derivative, respectively) diagram (see 
Fig. 1), most radio pulsars have a period of around one sec¬ 
ond and the period derivatives in the range of 10 to 10 
s s '. However, one group of sources having higher magnetic 
field (10*^“*^ G), relatively long periods (>1 s), and high pe¬ 
riod derivatives (10^*^ - 10“® s s“' ) are magnetars. The other 
group of sources are millisecond pulsars with old age (~ 10^ 
yr) and lower magnetic field (10* ® G). In the ATNF pulsar 
catalogue dHt], the longest spin period for normal radio pul¬ 
sars to date is ~8.51 sofPSRJ2144-3933 whose characteristic 
age is 2.72 x 10* yr and surface magnetic field is 2.08 x 10*^ 
G lH^. Except for short-lived glitches ll4^ . the rate of rota¬ 
tion of a NS drops steadily as a function of time, especially 
for millisecond pulsars which have very low values of P and 
very low timing noise. It is generally believed that the spin- 
down of the NS is due to magnetic dipole radiation (MDR) 
and electromagnetic torques in the magnetosphere M- 

Limits on DM from neutron stars have been deeply inves¬ 
tigated in a large amount of literatures so far 114^5Oil . Re¬ 
cently, Kouvaris and Perez-Garcia d have discovered that 
the electric charges could be expelled from the star as MCDM 
are accreted onto the NS. The escaping charged particles will 
provide an extra current, thus the torque is produced. They 
find that this mechanism yields an extra spin-down of neutron 
stars and the braking indices can be substantially smaller than 
3 predicted by the MDR model. 

We follow the philosophy of Kouvaris and Perez-Garcia 
iGtI], but explore how to constrain the s - m of MC particles 
via the currently observed pulsar periods. We consider the ef¬ 
fect of the additional torque by the extra currents due to the 
accretions of MCDM onto the NS on the period of pulsars. 
We also impose a limit of pulsar’s period in the P-P diagram 
on the period’s evolution equation as a cut-off. By doing so, 
the s - m of MC particles will be constrained under the given 
DM density. 

This paper is organized as follows. Firstly, we revisit the 
mechanism of the accretions of MCDM onto the NS (electro¬ 
magnetic accretion), then the rotational evolution equations 
of a NS considered the MDR and the torques, which are pro¬ 
duced by the Goldreich and Julian (GJ) current ifsil] /qj and 
the extra current /dm, are shown in Section 2. Secondly, we 
numerically show the bounds of the e - m of MC particles in 
Section 3. Finally, we present the conclusions in Section 4. 


II. THE MODELS 

First we show how the accretions of MCDM onto the NS 
can lead to an excess of currents ifT^ . If the Larmor radius 
n = of MCDM is much smaller than the curvature radius 
^curv = follow the magnetic field lines in the 

magnetosphere within the light cylinder of the NS and cross 
the surface of the star. v_i_ is the particle’s velocity perpendic¬ 
ular to the magnetic field B. As the accumulations of MCDM 
trapped in the NS, once the Coulomb repulsive force becomes 
larger than the gravitational force, the electric charges will be 
expelled from the star. The expelled charges are in the form 


of either electrons or protons, which depend on the charge of 
the accreted MCDM. Once the equilibrium between the ac¬ 
cretions of MCDM and the expulsions of electric charges has 
been established, stable extra currents will be formed. 

As stated previously, the slow-down of the rotation powered 
pulsars is due to the MDR and the braking torques provided 
by the outflowing plasmas if^ . The rotational kinetic energy 
loss rate of the NS is £ = IQ£l - -L, where I and Q are the 
moment of inertia and the angular velocity of the star, respec¬ 
tively, and L can be expressed as two components if^ 

£ — £orth sin 6 + £aiign cos 6. (1) 


In the above equation, £orth represents the “orthogonal” case 
that the angular momentum of the star is decreased by the 
MDR, Laiign represents the “aligned” case that the torque is 
produced by the electric current from escaping charged par¬ 
ticles that follow the open magnetic-field lines in the magne¬ 
tosphere of the NS, and 0 is the angle between the rotational 
and the magnetic axes (magnetic inclination angle). The two 
components can be written as 

“ 4 c 3 ’ “ 2c^ ’ 

where Bq and R are the magnetic field strength on the polar cap 
(the open magnetic-field lines in the light cylinder connect to 
the surface of the star) and radius of the star, respectively, c is 
the light speed, and fideath is the angular velocity below which 
the pulsar emission dies. The total current I can be written as 
the sum of Iqj and /dm, I - 7gj + ^dm, where Iqj represents the 
emission of relativistic charged particles from the surface cap 
regions (~ J^R^, Rc is the polar cap radius) of the NS, due to 
a particle acceleration gap potential Vgap of the order of 10'^ 
V if^ that develops along open magnetic field lines in the 
vicinity of the polar cap, and /dm represents the extra current 
from the expelled electric charges (electrons or protons) due 
to the accretions of MCDM onto the NS. 

According to the pioneering work of Goldreich and Julian 
iH, the GJ current is Iqj = nR^pcjc, where pcs is the GJ 
charge density. The polar cap radius can be expressed as Rc = 
/?sin6p = RiR/R^y^^, where Op is the angle of the polar cap 
to the center of the star and Rp = cPI2ji is the light cylinder 
radius Imf. In the polar cap region, the GJ charge density is 
given by l5l[] 


thus, we can easily deduce 

(4, 

For the extra current /dm due to the accretions of MCDM onto 
the NS, we directly adopt the results in Ref. 10, 


/dm ^ 1.0 X 10^‘^ee 


! Ppm GeV \/£\ 

\0.3GeVcm3/\ m /Is/ 


s-‘. (5) 
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FIG. 1: Distribution of all known pulsars on the P - P diagram. 
Dots are rotation powered pulsars (including normal pulsars 
and millisecond pulsars, which dominate the pulsar population) 
(from ATNF: http://www.atnf.csiro.aU/research/pulsar/psrcat/l. 
Triangles are magnetars (from McGill magnetar catalog: 
http://www.physics.mcgill.ca/~pulsar/magnetar/main.html I. The red 
star denotes the famous PSR J2144-3933, which has the longest 
spin period for normal radio pulsars. The dash lines of pulsar 
characteristic age are defined by PI(2P) and another dash lines 
of surface dipole magnetic field are conventionally defined as 
3.2 X W'^(PP)'^^, which are based on the MDR. 


Finally, using the Eqs. (1) and (2), we can get the rotational 
evolution equation of the NS 


dfl 

df 


4Ic^ 


sin^ 0 + 



^death \ 

Q. ) 



cos^ 6 


( 6 ) 

We define T = ^ for convenience. As we can see from Eq. 
(6), it easily recover the scene of the MDR for 4 = 0 (no ex¬ 
tra current). It is worth noticing that we can take «: 1 
safely, since normal pulsars couldn’t be older than 10^ to 10® 
years. We compare the evolution curve of period from Eq. 
(6) with that from the scene of pure MDR. Eig. 2 shows the 
differences between them. It is obvious from Eig. 2 that the 
differences become larger gradually as the increases of age. 
We use the 8.51 s period at 2.72 x 10® yr (which are the ob¬ 
servational data of PSR J2144-3933) as the cut-off of neutron 
stars, because it is the longest observed period of normal pul¬ 
sars and possible longer period should fall below the so-called 
“death line” on the P-P diagram, where pulse emission can’t 
be observed. As shown in Eig. 1, the location of PSR J2144- 
3933 on the P-P diagram indicates that the star could be a 
good candidate as the cut-off for the period evolution and the 
characteristic age of normal radio pulsars. However, both the 
millisecond pulsars and magnetars are unsuitable to constrain 
the e - m of MC particles. It is generally believed that the 
millisecond pulsars have a complicated and turned evolution. 
Recycling neutron stars to millisecond periods may be a key 
process. Additionally, the extra current /dm is much smaller 
than the corresponding GJ current /qj according to Eq. (4) for 


FIG. 2: The period evolution of a canonical NS with M = 1.4 Mq and 
R = 10 km for So = 10'^ G. The dash curve represents the spin-down 
induced by pure MDR (T = 0 in Eq. (6)). The solid curve shows an 
extra slow-down of the NS with magnetic inclination angle 6 = 45° 
by the additional torque due to the accretions of MCDM with m = \ 
GeV and s = 5.04 x 10“^. The labels for the red dots represent the 
8.51 s period at 2.72 x 10® yr. 


millisecond pulsars. The tiny differences have no imprint in 
Eig. 2 even if we overlook the recycled processes. Eor higher 
magnetic fields, there are only 28 currently known magne¬ 
tars and magnetar candidates 15^ . The physics of magnetars 
is unclear so far. Their emission is probably powered by the 
non-rotational kinetic energy or by their decay of super-strong 
magnetic field ifs^ . The pure MDR model may not be a 
good description for magnetars. Therefore, we think that it is 
seriously inaccurate and unreliable to constrain the parameters 
of MC particles by magnetars. 


III. NUMERICAL RESULTS AND DISCUSSIONS 

Eor convenience, we assume that the star consists of non¬ 
superfluid neutrons, protons and electrons, and suppose that 
the density of the star is uniform. In our calculations, we con¬ 
sider a canonical NS of M = 1.4 Mq and /? = 10 km, and take 
the initial period Pq - I ms, the moment of inertia of the star 
I - ^ 10"^^ g cm^, and the local DM density in the 

vicinity of the NS poM = 100 x 0.3 GeV/cm® (0.3 GeV/cm® 
represents the standard DM density around the Earth). 

The e-m phase space of MC particles could be constrained 
using the following strategy. Eirst, we plug a great number of 
groups of (e, m) into Eqs. (4), (5) and (6) for various magnetic 
field strength and different magnetic inclination angles under 
the given DM density in the vicinity of the NS. Then, we sim¬ 
ulate Eq. (6) up to 2.72 x 10® yr as the cut-off for the age of 
the NS. If the evolution of the rotational period of the star ex¬ 
ceeds the longest observed period of normal radio pulsars to 
date (~8.51 s), we consider this group of (e, m) invalid. 

The e-m phase space of MC particles constrained by a 











4 



S 


FIG. 3: The e - m phase space of MC particles constrained by a canonical NS of M = 1.4 Mq and I? = 10 km with Bo = 10*^ G for different 
magnetic inclination angles. 


canonical NS of M = 1.4 Mq and = 10 km with typi¬ 
cal Z?o = 10^^ G for different magnetic inclination angles is 
shown in Fig. 3. It can be seen from Fig. 3(a) that the charge 
of MC particles can reach up to e = 2.52 x 10“^ for heavier 
particles with m - \ GeV, the bound e > 2.52 X 10 ^ is ruled 
out apparently for MC particles with m-l MeV, and the limit 
e < 2.52x10“* is allowed for lighter particles with m - \ KeV. 
Furthermore, if the mass m decreases by one order of magni¬ 
tude, the charge e would reduce by the same one order of mag¬ 
nitude. As shown above, the lighter the mass, the smaller the 
charge. It is obvious from Eq. (6) that the rotational angular 
velocity of the star could decrease rapidly as A increases for 
the same magnetic field strength Bq, which result in the rise 
of the spin period even to the unreasonable value. Fig. 3(b), 
(c) and (d) show that the upper limit of s is inversely propor¬ 
tional to cos^ 6 for the same mass m of MC particles (see Eq. 
(6)). To summarize, we would present an upper limit of MC 
particles, ()xe<2.52xl 0“^/ cos^ 6, which is consistent 
with experimental and observational bounds ll^[^[33] . 

It is worth stressing that, as shown in Eq. (6), the rate of de¬ 
celeration of the rotation Q also depends on the magnetic field 
strength Bo- Eor Bo = lO" G, the bounds become weaker, the 
charge e of MC particles can’t be constrained for heavier par¬ 
ticles with m-l GeV, however up to s < 3.26 X 10 */ cos^ 0 
for m-l MeV and e < 3.26 x 10“®/ cos^ 0 for lighter parti¬ 
cles with m-l KeV. On the other hand, it can be seen from 
Eq. (5) that, the extra current /dm relies on the DM density 
Pdm in the vicinity of the NS. Although pdm around the NS is 
uncertain to date, we could draw a conclusion that, the larger 
the DM density, the smaller the charge for the MC particles 
with the same mass, e.g., if the DM density poM increases by 


one order of magnitude, the charge s would decrease by one 
order of magnitude correspondingly for the same mass m of 
MC particles. 


IV. CONCLUSIONS 

We have constructed the constraints of the s - m of MC 
particles by neutron stars based on an enhanced slow-down of 
neutron stars due to an extra current yield by the accretions of 
MCDM. Eor a canonical NS of M = 1.4 M© and R - 10 km 
with typical magnetic field strength Bq = 10*^ G, we have 
shown an upper limit of MC particles, (^-^f^) x e < 2.52 x 
10“^/ cos^ 0, which indicates the charge would be smaller as 
the mass of MC particles becomes lighter and is compati¬ 
ble with experimental and observational bounds. The specific 
limits are as follows, the charge of MC particles could rise to 
e - 2.52 X 10“^/ cos^ 0 for heavier particles with m-l GeV, 
the bounds e < 2.52 x 10“®/ cos^ 0 and e < 2.52 x 10“*/ cos^ 0 
are allowed for the MC particles with m-l MeV and m - 1 
KeV, respectively, i.e., if the mass m decreases by one order 
of magnitude, the charge e would reduce by the same one or¬ 
der of magnitude. However, for the NS with lower magnetic 
fields (10*®“" G), the bounds become weaker. In particular, 
for millisecond pulsars (10*“® G) and magnetars (10**“'® G), 
it is not appropriate to bound the e - m of MC particles. In 
addition, we have also investigated the influence on the e - m 
of MC particles for the different DM density pdm in the vicin¬ 
ity of the NS. It is obvious that, the larger the DM density, the 
smaller the charge for the MC particles with the same mass. 

The model we adopted is uniform stellar configuration. 
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However, it is well-known that the neutron stars constructed 
by the realistic equations of states can be approximated as the 
uniform case. Our work have shown an upper limit of the e-m 
of MC particles based on the new mechanism, and the results 
will be unchanged in the order of magnitude when considering 
the realistic equations of states. 
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